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ABSTRACT 

High-resolution simulations of cosmological structure formation indicate that dark 
matter substructure in dense environments, such as groups and clusters, may survive 
for a long time. These dark matter subhalos arc the likely hosts of galaxies. We examine 
the small-scale spatial clustering of subhalo major mergers at high redshift using 
high-resolution N-body simulations of cosmological volumes. Recently merged, massive 
subhalos exhibit enhanced clustering on scales ~ 100 — 300/i _1 kpc, relative to all 
subhalos of the same infall mass, for a short time after a major merger (< 500 Myr). 
The small-scale clustering enhancement is smaller for lower mass subhalos, which also 
show a deficit on scales just beyond the excess. Halos hosting recent subhalo mergers 
tend to have more subhalos; for massive subhalos the excess is stronger and it tends to 
increase for the most massive host halos. The subhalo merger fraction is independent 
of halo mass for the scales we probe. In terms of satellite and central subhalos, the 
merger increase in small-scale clustering for massive subhalos arises from recently 
merged massive central subhalos having an enhanced satellite population. Our mergers 
are defined via their parent infall mass ratios. Subhalos experiencing major mass gains 
also exhibit a small-scale clustering enhancement, but these correspond to two-body 
interactions leading to two final subhalos, rather than subhalo coalescence. 

Key words: cosmology:theory - methods:N-body simulations - galaxies: halos - 
galaxies : interactions 



1 INTRODUCTION 

A wealth of high redshift galaxy data is now accumu- 
lating, and many members in the resultant g alaxy zoo 
are t hought arise from galaxy mergers: q uasars llCarlberd 
[lJ39C|), Lyman Break Galaxies (LBG; see iGiavaliscol 20021. 
for r eview), submillimeter galaxies (SMG; see Blain et al.l 
20021, for review) , ultra- l umino us infrared galaxies (ULIRG; 
see ISanders fc Mirabell Il99fj . for review) and starburst 
and s tarburst remna nt galaxies (e.g., iBarnes fc Hernquistl 
Il99ll : | Noguchil |l99ll ) . Observational samples of these ob- 
jects at z > 1 are growing large enough to pro- 
duce statistical measurements of clustering which can be 
compared to candidates in numerica l simulations (e.g. , 
Giavalisco et al.lll998l; iBlain et alj|2004 ICroom et al.ll2005l: 
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Moore et al.ll 19991 ). These "subhalos" are thought to be the 



hosts of galaxies, and indeed the identification of galax- 
ies with subhalos reproduces many galaxy properties (e.g., 
Springel et all 1200 l |. l2005l: IZentner et "aT1l2005l : bower et al 



Ouchi et al.l 120051: ICoorav fc Ouchil 120061; iHermawi et al 



iKan 



2006: Kashikawa et al. 2006: iLee et al. 2006: IScott et al 



2006: lCoil et 3,1.1120071: iGawiser et al.l 20071: iS hcn ct al. 200, 
Francke et al.ll2008l : lMvers et alfcOOalYamauchi et allhoOE 
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In high resolution dark matter simulations, over- 
dense, self-bound, dark matter substructures in dense en- 



2006 ; IConrov et al.ll2006l : IVale fc Ostrikerll2006l : IWang et al 
20061 ). The complex dynamics of subhalos may thus be 
a good proxy for those of galaxies themselves, suggesting 
that galaxy mergers can be identified with subhalo merg- 
ers within simulations. We will use the terms galaxy and 
subhalo interchangeably hereon. Measurements of subhalo 
mergers can provide a quantitative reference for the iden- 
tification of merger related objects in observations and can 
aid in the correct interpretation of their clustering measure- 
ments. 

One question of particular interest is whether small- 
scale clustering can probe merger activity. There has been 
much recent discussion of a "small-scale clustering enhance- 
ment" of galaxy subpopulations in simulations and obser- 
vations, but several differing definitions of "enhancement" 
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exist, mostly stemming from different choices of reference. 
For instance, small-scale enhancement has been used to de- 
scribe clustering stronger than a power law extrapolated 
from larger scales, or clustering stronger than that of dark 
matter at small scales. However both of these behaviors 
are seen in non-merging samples. Luminous Red Galax- 
ies in SDSS are not thought to be associated with (re- 
cent) mergers, yet their correlation function is much steeper 
than the dark mat ter on scales of several hundred kpc 
l|Masiedi et al.ll2006l ). One expects objects which populate 
halos more massive than M* (the characteristic non-linear 
mass) will have an upturn in their correlation function on 
scales below the virial radius of the M* halos since at this 
scale the clustering is dominated by pairs of objects within 
the same h a lo and the mass functi on is very steep (e.g., 
ISeliakl 1200(3 : iPeacock fc Smith! 120001 ). Quasars (thought to 
be associated with mergers) do exhibit a sharp upturn in 
clustering on very small scales (25 — 50 h ~ kpc) at z ~ 1 — 3 
l|Hennawi et al.ll2006l ; iMvers et al.ll200Sl ), but a comparison 
with galaxy clustering measurements at these scales and red- 
shifts is lackingQOn slightly larger scales, low z quasar clus- 
tering observations show n o excess above a power law (e.g. 
IPadmanabhan et al.ll200"il . and discussion therein). 

On the other hand, merger related clustering effects 
are not unexpected. Recently, various authors have shown 
that the large-scale clustering of dark matter halos depends 
on their formation histories (known as "assembly bias"; 



Sheth fc Tormenll2004l ; lGao et al.H 2005: Wcc hsler et al l 



12006 ) , and in particular, on recent halo merger activity 
("me r ger b ias" ; IScannapieco fc Thackerll2003l ; IWetzel et al.l 



( me r ger b ias' 
12001 [2008)0 



Any such history or merger dependent clus- 
tering breaks the usually assumed direct link between 
large-scale clustering amplitude and halo mass. Analyti- 
cal modeling of the small-scale clustering of quasars has 
been compared to observations, assuming that quasars are 
mergers and that m ergers occur in denser environments 
jHopkins et al.ll2008T l. 

In this paper, we focus on the small-scale clustering of 
subhalo mergers and its relation to their host halos. We use 
dark matter simulations to compare the small-scale clus- 
tering of recently merged subhalos at high redshifts to the 
clustering of the general subhalo population of the same 
infall mass. To interpret our results, we relate the merged 
subhalos to their host dark matter halos using the formal- 
ism of the halo model (ISeliakll2000l ; IPeacock fc Smithll2000l ; 
ICoorav fc Shethll2002l ). examining the dependence of sub- 
halo mergers both on their host halo masses and their halo 
radial distribution profiles. 

In fj2] we briefly describe the simulations and summa- 
rize the main properti es of subhalo trackin g and merger def- 
initions as detailed in IWetzel et all ([2009), hereafter called 
Paper I. In Sj3]we compare the small-scale clustering of the 
recently merged and full population of subhalos. We also de- 
scribe results for the alternate mass gain merger definition. 
In !j4] we use the framework of the halo model to identify 



1 It is not clear how to interpret a galaxy clustering measurement 
on scales s maller than the galactic radiu s. 

2 Sec also Furlanctto & Kamionkowski (2006) for analytic esti- 
mates. [Perciva l et al.l j2003T) for a simula tion limit on the effect, 
ICroton et al.l d2007h ; iTinker et al.l fcOOgT) for halo assembly bias 
applied to galaxy clustering. 



contributions to merger clustering and to understand the 
relation of subhalo mergers to their host dark matter halos. 
We note properties of merger pairs in iJS] and we summarize 
and discuss our results in 36] 



2 NUMERICAL TECHNIQUES AND MERGER 
DEFINITIONS 

2.1 Simulations and Subhalo Tracking 

Our simulation and subhalo finding and tracking details 
are discussed extensively in Paper I; we summarize the 
main feat ures here. W e use two dark matter only N-body 
TreePM (|Whitell2002T ) simulations of 800 3 and 1024 3 par- 
ticles in periodic cubes with side lengths 100/i _1 Mpc and 
250 /i _1 Mpc, respectively. For our ACDM cosmology (£7 m = 
0.25, = 0.75, h = 0.72, n = 0.97 and < t 8 = 0-8), in agree- 
ment with a wide ar r ay of observations (ISmoot et al. 19921; 
i Tegmark et al.ll2006l ; iReichardt et al.ll200Sl ; iKomatsu et all 
2009), this results in particle masses of 1.4 x W s h- 1 M Q 
(1.1 x 10 9 '/i -1 Mq) and a Plummcr equivalent smoothing of 
4 ft -1 kpc (9 ft -1 kpc) for the smaller (larger) simulation. Out- 
puts were spaced every 50Myr (~ 100 Myr) for the smaller 
(larger) simulation, from z ~ 5 to 2.5. Additional outputs 
from the smaller simulation were retained at lower redshift, 
spaced every ~ 200 Myr down to z = 0.6, below which we 
no longer fairly sample a cosmological volume. 

We find subhalos (and sometimes subhalos in subhalos) 
by first generating a cata log of halos using the Friends-of- 
Friends (FoF) algorithm (iDavis et al.l lT985) with a linking 
length of b = 0.168 times the mean inter-particle spacing[j 
We keep all groups that have more than 32 particles, and 
halo masses quoted below are these FoF masses. Within 
these "(host) halos" we then identify "subhalos" as grav- 
itationally self-bound aggregations of at least 20 particles 
bounded by a density saddle point, using a new imple - 
mentation of the Sub find algorithm (|Springel et al.l [200lh ■ 
The central subhalo is defined as the most massive subhalo 
within its host halo, and it includes all halo matter not as- 
signed to satellite subhalosQ Thus, the central subhalo mass 
represents halo mass not bound to any satellite subhalos and 
so it should track the gas mass available to accrete onto the 
central galaxy. Subhalo positions are those of their densest 
particle, and the halo position is that of its central suhalo. 

Each subhalo is given a unique child at a later time, 
based on its 20 most bound particles. We track subhalo his- 
tories across four consecutive outputs at a time since subha- 
los can briefly disappear during close passage with another 
subhalo, i.e. "fly-by's". Additionally, at these redshifts, the 
distinction between a central and satellite subhalo is often 
not clear-cut, especially for halos undergoing rapid merger 
activity which are highly disturbed and aspherical. In partic- 
ular, our tracking also can produce "switches" : if a satellite 



3 The longer linking length of b = 0.2 is often used, but it is more 
susceptible to joining together distinct, unbound structures and 
assigning a halo that transiently passes by another as a subhalo. 
Halos at our mass and redshift regime have roughly a ~ 15% 
lower FoF mass than for b = 0.2. 

4 Since SubBnd defines the central subhalo as the most massive 
subhalo, it is not necessarily the case that the central subhalo 
contains the most dense or deepest potential particles. 
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becomes more massive than the central, it becomes the cen- 
tral while the central becomes a satellite, often switching 
back in the next output Q. In these cases, there is typically 
not a well-defined single peak that represents the center of 
the halo profile. 

We assign to a subhalo its mass when it fell into its 
current host halo, i.e., its subhalo mass when it was last a 
central subhalo, Mi n f. Centrals are assigned Mj n f as their 
current bound mass. Subhalo infall mass has been shown 
to correlate with galaxy stellar mas s i|Vale fc Ostrikerll2006l ; 
IWang et al]|200rj ; (Yang et all 120091 )1^1 If two satellite subha- 
los merge, their satellite child is given the sum of their infall 
masses. During a switch, when a satellite becomes a central, 
it would acquire all halo mass not bound in other subhalos. 
To avoid strong fluctuations in Mi n f, a central subhalo is 
assigned its halo's current self bound mass only if it was a 
central in the same halo in the previous output. If a central 
was a satellite or a central in a different (smaller) halo in the 
previous output, it is assigned the sum of its parents' Mi n f. 

We select subhalos with M inf > lO 12 ft _1 M in the 
larger simulation and scale down to M ln t > 10 11 h -1 Mq 
in the smaller, higher resolution simulation, by requiring 
consistency between the two simulations in their overlap 
regime^ Halos of mass cross below M t , 

the characteristic mass of collapse, at z = 1.5 (z = 0.8), 
so we probe M > M* subhalos for much of the red- 
shift range we consider. We also expect our sample of 
Mint > 10 12 h~ x Mq subhalos to approximately correspond 
to L > L« galaxies at the r edshifts we examine (see, 
e.g., IConrov fc Wechslerl (I2008T I for halo-galaxy mass re- 
lation based on abundance matching, but note that they 
match stellar mass to halo mass, which is typically 10— 15% 
higher than subhalo infall mass) . Additi onally, most massive 
galaxies are gas-r i ch (blue) at z > 1 (|Cooper et al.l 120071 : 
iGerke et al. I l2007l ; iHopkins et ai]|2008l ). possessing enough 
gas to be actively star forming. Thus, we anticipate that 
most, if not all, mergers we track have the capacity to drive 
galaxy activity such as starbursts and quasars. 



2.2 Merger Criteria 

We select a subhalo as a major merger (henceforth merger) if 
its two most massive parents, with Mi n f,2 < Mjnf.r, satisfy 
Minf,2/M"inf,i > 1/3. As mentioned above, galaxy mergers 
with stellar mass ratios closer than 3:1 are expected to drive 
interesting activity, e.g. quasars and starbursts. Unless oth- 
erwise stated, we use the shortest simulation output spacing 
to define the merger time interval, corresponding to 50 Myr 



5 In more detail, a switch occurs when the density peak of a satel- 
lite (above the background) contains more mass than is within the 
central subhalo's radius at the position of the satellite 

6 As has subhalo maximum circular velocity at infall, V c i n f 
llConrov et al.ll2006l; iBerrier et alj|2006h . We compare M inf and 
Vq inf i n detail in Paper I. 

7 Using the Mill e nnium simulation jSpringel et al.l |2005|) . 
iKitzbichler fc White! 11200811 require an analytic model for satellite 
infall times after subhalo disruption to match small-scale galaxy 
clustering at z ~ 0. We do not expect this numerical disruption 
to significantly bias our results since our 100 /i — x Mpc simulation 
has higher mass and temporal resolution. 



(~ 100 Myr) for M inf > 10 11 (1O 12 )/i _1 M at z > 2.5, and 
~ 200 Myr for all masses at z < 1.60 

Other definitions of mergers produce significantly dif- 
ferent merger samples. In re l ated work, which inspired our 
investigation, iThacker et alj l|2006h used a dark matter plus 
hydrodynamic simulation to measure the small-scale clus- 
tering of subhalos with recent large mass gains, finding that 
these subhalos have enhanced small-scale clustering relative 
to a population with the same large-scale (> l/i _1 Mpc) 
clustering. Mass gain is convenient in that it does not re- 
quire histories beyond the previous time step, and mass gain 
is unambiguously defined for all subhalos. However, using 
our simulations and subhalo finder, the resulting sample is 
almost entirely different from the one defined above0 Specif- 
ically, using a mass gain merger definition in our simulations 
led to 'mergers' where the two contributing galaxies almost 
always remained as distinct entities after the merger event. 
The most common instance of major mass gain is a satellite 
subhalo gaining mass during its movement within its host 
halo, particularly as it m oves away from the h alo center 
(see Fig. 2 in Paper I and lDiemand et al.1 ((2007) for exam- 
ples). A subhalo can also gain mass by stripping material 
from the outskirts of a nearby subhalo. In 75% of the cases 
of major mass gain, one of the progenitors contributed less 
than 10% of its mass to the resulting 'merged' child. The 
most bound particles (where we expect the stellar compo- 
nent of a galaxy) were unaffected. We did find significantly 
increased small-scale clustering for these mass gain subha- 
los, attributable to the remaining nearby subhalo which just 
'merged' with it. Simil ar issues in using mass g ain to define 
mergers were noted in lMaulbctsc h~et al.l (|2007T l. 



3 SMALL-SCALE SPATIAL CLUSTERING 

When examining the effects of recent mergers on spatial 
clustering it is important to define an appropriate compar- 
ison sample. We have chosen all subhalos above the same 
given infall mass cut as the mergers, with a matched Minf 
distribution. We match Minf to remove any possible artificial 
biasing from, e.g., mergers preferentially occurring for sub- 
halos of higher mass. Using only a mass cut without match- 
ing the mass distribution leads to a similar, but weaker, 
effect. If infall mass is a good proxy for stellar mass, our 
merger and comparison samples correspond to populations 
matched in stellar mass. 

To measure the small-scale (~ 100 — 1000 /i kpc) 
clustering of subhalo mergers relative to the clustering 
of the general subhalo population, we measure the cross- 
correlation function^ of the merged and general population, 
£ mg , and the auto-correlation function of the general popu- 
lation, £ 99 . Our limited volume unfortunately does not allow 



8 This gives 280 (490) mergers at z = 2.6 (z = 1) above the 
lower mass cut in our smaller, higher resolution simulation. 

9 For M inf > W 11 h~ 1 M (3 subhalos at Z = 2.6, 492 have 
-Mcont,2/Af con t,l > g> compared to the 260 for our infall mass 
ratio definition, with only 7 overlapping the two sets. 

10 Not only is the cross-correlation of the mergers with the gen- 
eral population interesting in itself, but it also provides better 
statistics. 



4 Wetzel, Cohn & White 




0.1 1 10 0.1 1 10 

r [h-'Mpc] r [h-'Mpe] 

Figure 1. Cross-correlation of recently merged subhalos with all subhalos, § ms (r), and the auto-correlation of all subhalos, ^ gg {r), at 
z = 2.6, for M; n f > 10 12 h~ 1 Mq (left) and M; n f > 10 11 h —1 Mq (right). Infall mass, M; n f, is matched between the merged and full subhalo 
samples. Top Left: £ mff (r) for subhalos merging within the last 130 Myr (long-dashed) and !; gg (r) for all subhalos (solid). Bottom Left: 
Ratio of the cross- and auto-correlations above, for mergers within the last 130 Myr (long-dashed) and 500 Myr (short-dashed). Top 
Right: 5 m9 (r) for subhalos merging within the last 50 Myr (long-dashed) and £, gg (r) for all subhalos (solid). Bottom Right: Ratio 
of the cross- and auto-correlations above, for mergers within the last 50 Myr (long-dashed) and 250 Myr (short-dashed). Higher mass 
subhalos show stronger enhanced clustering from mergers, but at both masses no signal persists for subhalos > 500 Myr after merging. 
Errors are given by yj N ps .i v and do not include sample variance. 



us to sub-divide our simulation to measure the sample vari- 
ance error. We show yj N ps ,i T errors on the correlation func- 
tion points, but we caution that this may underestimate the 
error by up to a factor of 2. Our clustering measurements 
are limited on small scales by the force resolution and on 
large scales by the simulation volume. We present results on 
scales where these effects are minor, and we further discuss 
the effects of finite simulation volumes and the statistics of 
massive halos in H4.3I 



Figure [T] shows the spatial clustering of recently merged 
subhalos and the general subhalo population at z = 2.6, for 
two mass regimes and several merger time intervals. These 
results are representative of our results at other redshifts. 
Over the smallest time intervals, both high and low mass 
subhalos have a rise and a decline relative to the general 
population, with the rise being most prominent for higher 
mass subhalos and the decline most prominent for lower 
mass subhalos. For Af Jnf > 10 11 ft _1 M , we find an up- 
per limit of 1.8x enhancement at 70/i _1 kpc, increasing to 
3x that of the general population at ~ 150 h~ kpc for 
Mjnf > 10 12 hT 1 Mq. In addition, lower mass subhalos ex- 
hibit a deficit at 100 — 300 /i _1 kpc. The enhancement/deficit 
declines rapidly with the time since the merger, and we see 
no signal > 500 Myr after the merger. As time progresses 
central mergers will become satellites in larger halos and 
satellite mergers will move within their host halos (and per- 
haps merge with the central), washing out the correlation 
between the merger and its halo properties that we see at 
the time of the merger. 



4 HALO OCCUPATION DISTRIBUTION AND 
RADIAL PROFILE 

The halo model dPeacock fc Smith! l200d : ISeliakl l200d : 
ICoorav fc Shethll2002n " can provide insight into the observed 
clustering signals of the merged subhalo population in Fig. [T] 
In this framework, galaxies populate dark matter halos such 
that their large-scale spatial clustering is determined pri- 
marily by the clustering of their host dark matter halos ( "2- 
halo term"), while their small-scale spatial clustering arises 
from galaxies in the same host halo ("1-halo term"). The 
objects occupy dark matter halos according to a Halo Oc- 
cupation Distribution (HOD) and have some radial profile 
within these halos. Since the clustering of recently merged 
galaxies differs from that of all galaxies, we expect mergers 
to differ from the general galaxy population in their HOD 
and/or profile. 



4.1 HOD of Subhalos 

Figure [2] shows the HOD of subhalos (central and satellite) 
at 2 = 2.6 (corresponding to the populations in Fig. [1} for 
all subhalos and for all subhalos within a halo containing 
a recently merged subhalo. We see clearly that halos host- 
ing subhalo mergers tend to have more subhalos. This in- 
crease occurs for both our high and low mass samples and 
for all redshifts we probe. For massive subhalo mergers the 
increase is larger and tends to rise to larger halo mass, which 
enhances the cross-correlation for mergers. Although lower 
mass subhalo mergers also have more subhalos per halo, the 
relative number does not change strongly with increasing 
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Figure 2. Top Panels: Halo Occupation Distribution (HOD) at z = 2.6 for all subhalos (solid) and subhalos in a halo hosting 
a recently merged subhalo (dotted). Shown are subhalos with Af; n f > l0 12 h~ 1 M(? > and mergers occuring within 130 Myr (left), and 
M; n f > lO 11 ^ - 1 Mq and mergers occurring within 50 Myr (right). Bottom Panels: Ratio of the above HOD's of subhalos in a halo 
hosting a merger to that for all subhalos. There are more subhalos in halos with recently merged subhalos, and the effect is stronger 
for more massive subhalos. For recently merged massive subhalos, which show a strong increase in small-scale clustering, the relative 
subhalo excess in halos with mergers tends to increase with halo mass. Similar trends persist at all redshifts we probe. 



halo mass - host halo mass does not significantly influence 
merger statistics. 

Note that we expect some increase in the number of 
subhalos above a given mass threshold from mergers of sub- 
halos just below the mass threshold, while mergers between 
subhalos above the threshold decrease the number of subha- 
los. Which effect wins out requires detailed simulations such 
as ours. 



4.2 Central and Satellite Cross-Correlation 

We now distinguish between contributions from satellite and 
central subhalos to the correlation function and HOD of 
mergers and the general population. There are inherent sub- 
tleties in this breakdown; as mentioned earlier, the identifi- 
cation of satellite vs. central subhalos is not entirely clear- 
cut at these masses and redshifts. In particular, one type 
can switch to another, and often does for mergers (see Pa- 
per I for more detail). Also, other tracking schemes might 
alter the relation between central and satellite assignments 
we use. If switched satellite mergers are assigned as centrals, 
this will of course not change the observed clustering signal 
but would alter where the contribution from our switches 
shows up in the breakdown of central and satellite effects. 

Figure [3] shows the ratios of the cross-correlations for 
central mergers and satellite mergers to the auto-correlation 
of the full sample (the left hand side of Fig. [T] above), for 
M in f > 1O 12 /i _1 M subhalos at z = 2.6. The ratios of the 
cross-correlations of the full central and satellite populations 
to the full subhalo population are also shown. As before, 
the full population is matched in M[ n { to the merged pop- 
ulation. In general, satellites have a larger cross-correlation 
with all subhalos compared to centrals because the mini- 
mum host halo mass for a subhalo of a given mass is larger 



M,„ t >lel2 h-' M, 



all satellites 

merged satellites 

all merged 

merged centrals 

all centrals 




r [h-' Mpc] 

Figure 3. The cross-correlation of the component popu- 
lations (merged and all satellites, merged and all centrals) 
with all subhalos (the latter with M; n f matched distribu- 
tion), at z = 2.6 for A/ inf > 10 12 /i _1 Mq, correspond- 
ing to the left hand panel of Fig. [T] Top to bottom at 

far left are ^sat-ail/^aili ^sat.mcrgcd-all/Calli ^all.mcrgod-all/^all , 
Seen, merge d-all/Call: sall-all/^all i 5con- all/^all • 



for a satellite than a central subhalo (increasing the large- 
scale clustering) , and a satellite will always have a central in 
its halo but not vice versa (increasing the small-scale clus- 
tering). The satellite and central cross-correlations must be 
summed, weighted by their population's fraction of the full 
population to get the full cross-correlation. 

The merged central and merged satellite contributions 
to the cross-correlation differ from their counterparts in the 
full population. Generally, the merged satellites have de- 
creased small-scale clustering, relative to all satellites, while 
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merged centrals have enhancement relative to all centrals. In 
addition, satellites comprise a larger fraction of the merger 
population than of the full subhalo population (typically 
~ 1.5-2X larger, see Paper I), i.e., A r S at,mo rgc d/A r con,mergcd / 
A?sat /-/Veen , changing the relative weights satellite and cen- 
tral contributions in the full cross-correlation[3 These three 
trends persist at all redshifts and subhalo masses, but their 
relative strengths vary to give the different behavior seen in 
Fig. Q] for different subhalo mass cuts. 

4.3 Central and Satellite HOD 

We can employ this central/satellite split for the HODs as 
well. Figure U shows the HODs for two mass ranges and two 
redshifts: at left is the HOD for high-mass subhalo merg- 
ers at z = 2.6, corresponding to Fig. [T] (left), and at right 
is the HOD corresponding to the less massive subhalos at 
z — 1, for a longer (230 Myr) merger time interval. The 
cross/auto-correlation ratio for this latter sample is similar 
to that shown in Fig. [T] (right) , but we show the results at 
lower redshift to illustrate the trends we see over a larger 
host halo mass range. 

Figure [4] (left) shows that high-mass recently merged 
central subhalos occupy halos with an excess of satellites 
(relative to the average) at most host halo masses (dot- 
ted curve). The physical extent of the enhancement, ~ 
300/i~ 1 kpc, coincides with the virial radius of the most 
massive halos (~ 4 x 10 13 /i _1 M Q ) in our simulations at 
z ~ 2.6. The enhancement of the satellite population in ha- 
los hosting recent central mergers is not as large for lower 
mass subhalos (the curve for this quantity is almost indistin- 
guishable from that of the full satellite population, thus it 
is not shown in Fig. fright). This corresponds to the weak- 
ened enhancement, and decrement, in Fig. [1] In terms of 
the central/satellite subhalo breakdown, the increased satel- 
lite number for recently merged central subhalos has the 
strongest correspondence with increased small-scale cluster- 
ing. 

The other trends that we see appear across our redshift 
and mass regimes. For a fixed subhalo mass cut, the fraction 
of centrals that have experienced a recent merger increases 
with halo mass. This causes an enhancement of the central 
merger cross-correlation with the full population, since the 
rise in central merger fraction with halo mass gives the over- 
all recently merged central population proportionally more 
satellites when all halo masses are summed. 

It is worth noting two points about the relative increase 
of central mergers in higher mass halos. The increase itself 
might be surprising, as both the number of satellites and the 
dynamical friction timescale increase approximately linearly 
with halo mass. One expects that the increase of possible 
satellites and the slowdown of their approach to the center 
would then cancel, rather than producing increased num- 
bers of central mergers. Looking instead at merger parent 
types, we found that the central-satellite merger HOD does 
not increase as steeply as the merged central HOD. Cen- 
tral mergers in the highest mass halos are often satellite- 

11 Since satellites preferentially reside in higher mass halos, if 
the recently merged population is comprised of a higher fraction 
of satellites, it will also exhibit enhanced large-scale clustering. 



satellite mergers whose child is a central (switches), rather 
than central-satellite mergers. While this is a small fraction 
(3%) of all central mergers, it is a larger fraction of those in 
high-mass halos. 

Secondly, the increase of central mergers to higher mass 
halos suggests a possible error from neglecting high mass 
halos which are too rare to occur in the simulation volume. 
Their effect can be estimated analytically by extrapolating 
the central merger fraction as a function of halo mass, using 
the halo model to find the relative contribution of merged 
and all centrals to the cross-correlation at some given radius, 
and seeing how this changes as larger masses are included. 
We find that, on average, neglecting halos more massive than 
those found in our simulations only causes a small change 
in the cross-correlation of centrals with the full population. 
However, our largest halos are quite rare, and so the merger 
fraction in them does not always tend to the average. The 
rise in central merger fraction with halo mass is largest for 
our Mi n { > 10 11 h -1 Mq subhalos and thus of most concern. 
As a check, we averaged our cross-correlation for the smaller 
simulation over several outputs to confirm the trends we see 
in Fig. [T] (right). The output to output variation does average 
to this trend, but individual outputs can show different sign 
effects, albeit with significant error bars. This is especially 
true at lower redshift, where our imperfect sampling of high 
mass halos plays a larger role. It might be counterintuitive 
that 300/i -1 kpc clustering is not very well measured in a 
100 /i _1 Mpc box, but the peak of the power spectrum is at 
large scales in a ACDM model and the scatter in the cross- 
correlation is driven by the contribution from rare, massive 
halos which host many satellites. 

We now turn to the satellite HOD. Figure [4] (bottom) 
shows that the fraction of merged satellites to all satel- 
lites decreases with increasing host halo mass (long-dashed 
curve). However, a significant fraction of merged satellites 
are satellite-central switches, which dominate in low-mass 
halos where a central and satellite are more likely to be com- 
parable in mass and thus able to switch. Additionally, some 
merged satellites are central-satellite mergers which then fell 
into the host halo. Examining instead merger parent type, 
satellite-satellite mergers are essentially a constant fraction 
of satellite subhalos across all halo masses (short-dashed 
line). We see a slight increase for massive host halos at 
2 = 2.6 and no increase at z — 1. Additionally, we see weak 
evidence for a cutoff in the satellite-satellite merger HOD 
for the most massive halos, consistent with the idea that 
increasing relative velocities of satellites in the most mas- 
sive h alos cuts off satellite-satellite mergers jMakino fc HutJ 
I1997I V However, the high-mass halo statistics are poor in our 
modest simulation volumes. 

A simple argument can be made which reproduces the 
observed scaling of the number of satellite-satellite merg- 
ers with the overall satellite population, based upon the 
HOD scaling with increasing halo mass/volume. For satel- 
lites with random orbital parameters in a halo of a fixed 
mass, the satellite-satellite merger probability will scale as 
number density squared (n^t), so the fraction of satellites 
undergoing satellite-satellite mergers scales as n sat . For ha- 
los above a few times the mass of a given satellite pop- 
ulatio n, the number of sa tellites scales with the halo mass 
(e.g- jKravtsov et al.ll2004h . which scales approximately with 
the halo volume, so AT sat tx Af ha io oc 14aio- Thus, n sat « 
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Figure 4. Halo Occupation Distribution (HOD), broken down in terms of merger types: centrals, satellites, recently merged centrals, 
and recently merged satellites. Top Left: Subhalos with Af; n f > 10 12 h~ 1 Mq at z = 2.6, corresponding to the left hand side of Fig. [T] 
and the component contributions in Fig. \3\ The merger time interval is 130 Myr. Also shown is the HOD for satellites in halos with a 
recently merged central. Bottom Left: Ratio of HOD's of recently merged centrals to all centrals and of recently merged satellites to 
all satellites. Top Right: Subhalos with M in f > lO 11 /! - 1 Mq at z = 1 and mergers occurring within 230 Myr, with cross-correlation 
similar to the 250 Myr interval in Fig. [T] right at z = 2.6. The HOD is also shown for satellite-satellite mergers. No strong excess is seen 
for satellites in halos with a recently merged central (not shown). Bottom Right: Ratio of HOD's of recently merged centrals to all 
centrals, recently merged satellites to all satellites, and satellite-satellite mergers to all satellites. Also shown is the ratio of HOD's of 
subhalo mergers to all subhalos (regardless of type). Satellite-satellite mergers, and subhalo mergers (regardless of type), show no strong 
dependence on halo mass, though our statistics are not sufficient to rule out a cutoff at the highest mass. 



constant, so the satellite-satellite merger fraction remains 
roughly constant. 

Finally, the lower right panel of Fig. U shows that when 
we do not split by type, the subhalo merger fraction is essen- 
tially independent of halo mass. Again, this contradicts the 
frequently expressed intuition that mergers are less frequent 
in higher mass halos, though we caution that our statistics 
are poor for cluster-mass halos, and we are working at high 
redshift. 



4.4 Radial Distribution Profile 

Figure [5] shows the other ingredient required to predict clus- 
tering: the (stacked) satellite radial distribution profile, for 
subhalos with Mj nt > 10 11 }i~~ 1 Mq at z = 2.6 and a merger 
time interval of 50 Myr. We choose this mass regime and 
redshift because the high temporal resolution allows us to 
accurately identify the locations of mergers, though the be- 
havior across all our mass and redshift regimes is consis- 
tent with these results. Shown are the profiles for all satel- 
lites, recently merged satellites, recently merged satellites 
from satel lite-satellite parents (because of switches), and 
the NFW (|Navarro et al.ll 19961) halo density profile. We nor- 
malize the density distributions to the average value within 
We use 2r2ooc because an appreciable fraction of 



1 We calculate the halo virial radius, r200o i- e - the radius 
within which the average density is 200 X the critical density, from 
the FoF (b = 0.168) mass by first converting to M20OC assum- 



satellite mergers occur just outside r2ooc given the highly as- 
pherical geometry of halos at this mass and redshift regime. 
The satellite subhalo profile well-traces the NFW profile 
from r2ooc down to 0.2r2ooc, and shows no dependence on 
subhalo mass, consistent with the dist ribution of subha- 
los se lected on infall mass at z = of iNagai fc Kravtsovl 
(2005)0 The deficit below 0.2r2ooc is a result of finite spa- 
tial resolution, tidal disruption of subhalos near the halo cen- 
ter, and subhalos temporarily disappearing from the sample 
during a fly- by (see Fig. 1 of Paper I). 

As Fig. [5] shows, the radial distribution profile of satel- 
lite mergers approximately follows that of the entire satel- 
lite population. The profile can be thought of as a measure 
of how satellite mergers correlate with centrals, bearing in 
mind that the distances are scaled by the halo virial radius. 
Relative to the density distribution of all satellites, recently 
merged satellites have a slightly more concentrated profile, 
with enhanced probability of being at 20% — 40% of the 
virial radius (long-dashed). This produces a relative deficit 
in clustering of recently merged satellites out to the halo 
virial radius. Considering instead recently merged centrals, 
we find that the profile of satellites in halos with a recently 



ing a spherical NFW density profile, and then taking M20OC = 
200f p c r 2 3 00c . 

13 When selecting subhalos on instantaneous bound mass, more 
massive halos are preferentially biased to the outer regions of a 
halo, since mass strip ping leaves few m assive subhalos in the inner 
region of a halo fsee lNagai fc K ravtsov 20051). 
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Figure 5. Top: Halo radial distribution profile of satellites with 
M inf > 10 11 h^M® at z = 2.6, for all satellites (solid), recent 
mergers (long-dashed), and mergers from satellite-satellite par- 
ents (short-dashed). The merger time interval is 50Myr. Also 
shown is the halo NFW density profile (dot-dashed), using an 
NFW concentration of c = 4, typical for halos at this mass and 
redshift. Density distributions are normalized to the average den- 
sity of the population within 2r200c- Bottom: Ratio of the re- 
cently merged satellite normalized density to that of all satel- 
lites (long-dashed) and the same for recent mergers from satellite- 
satellite parents (short-dashed). We find no dependence of these 
results on satellite infall mass or redshift. 



merged central traces the general satellite population with- 
out significant deviation (not shown). 

While recently merged satellites exhibit a more concen- 
trated profile, mergers between satellite-satellite subhalos 
preferentially occur in halo outskirts and are comparatively 
less common in the central regions (short-dashed curve). 
Thus, the enhanced probability of finding recent satellite 
mergers at small scales is driven by switches (when a satel- 
lite merges with a central and results in a satellite). 



5 MERGER PAIRS 

While we have focused on measuring merger clustering via 
the cross-correlation of mergers with the general popula- 
tion, we can also measure directly the statistics of our lim- 
ited number of merger pairs and distinguish the different 
origins for the pair members. We focus on pairs of subha- 
los within 250 /i _1 kpc which have both undergone a merger 
within the last ~ 250 Myr. In all such cases of close merger 
pairs, both subhalos inhabit the same halo. These pairs are 
quite rare: at z = 2.6, there are 24 (2) per (100 /i _1 Mpc) 3 
for Mm > 10 11 (10 12 )/i" 1 M Q , while at z = 1 there are 2 



merger pairs are satellite-satellite subhalos. In all cases, the 
central-satellite pairs arise when a satellite-central merger 
occurs simultaneously with a satellite-satellite merger in a 
single halo. For the rarer cases of satellite-satellite merger 
pairs, two-thirds of the recently merged satellites arise from 
satellite-satellite parents within a halo, and one third arise 
from a satellite-central merger (a switch). 



6 SUMMARY AND DISCUSSION 

Using high-resolution dark matter simulations in cosmolog- 
ical volumes, we have measured the small-scale spatial clus- 
tering for massive subhalo mergers at high redshift and com- 
pared against the clustering of the general subhalo popula- 
tion of the same mass. We have described the merger pop- 
ulations in terms of their HOD and radial profile, including 
the contributions of centrals vs. satellites. We assign subha- 
los their mass at infall, as a proxy for galaxy stellar mass, 
but make no further attempt to model the baryonic com- 
ponent. We consider mergers with < 3:1 infall mass ratios, 
motivated by the expectation that these can trigger activity 
such as quasars or starbursts. Our main results are: 

• At z = 2.6, recently merged, massive (M- m f > 
10 12 h~ 1 MQ) subhalos exhibit enhanced small-scale cluster- 
ing compared to random subhalos with the same infall mass 
distribution. This excess peaks at 100 — 300/i _1 kpc, while 
the clustering exhibits a dip at slightly larger scales. Lower 
mass subhalos (Mi n f > 10 11 K~ 1 Mq) exhibit signs of a small 
rise in clustering at < 100/t -1 kpc, with a deficit at slightly 
larger scales, though our results are noisy in this regime. We 
find similar behavior at z = 1. The merger signal weakens 
rapidly with time, vanishing for time intervals longer than 
500 Myr after the merger. 

• Considering their HOD, halos hosting recently merged 
subhalos tend to have more subhalos. This enhancement is 
stronger for more massive subhalos and exhibited growth 
with halo mass. 

• Breaking the contributions to the cross-correlation into 
those from satellite and central subhalos, the recently 
merged massive centrals, which show the largest enhance- 
ment of small-scale clustering, preferentially occupy halos 
with more satellites. More generally, recently merged cen- 
trals occupy higher mass halos and recently merged satellites 
occupy lower mass halos with fewer (perhaps no other) satel- 
lites. The resulting increase for the central cross-correlation, 
decrease for the satellite cross-correlation, and change in the 
ratio of their contributions combine to produce the cross 
clustering enhancement or decrement. 

• For the range of halo masses we probe, the (type- 
independent) subhalo merger fraction is independent of host 
halo mass. We find similar behavior for the halo mass depen- 
dence of satellite-satellite mergers, i.e. the fraction of satel- 
lites that experience a merger with another satellite does not 
depend on halo mass, except perhaps in our poorly sampled 



per (100 /i _1 Mpc) 3 for M inf > 10 11 /i _1 M a and none at the largest mass halos (~ 10 14 /i _1 M Q ) 
higher mass. 

For two subhalos within the same halo, the two types of 
possible pairings are satellite-central and satellite-satellite. 
For M inf > 10 12 ft -1 Afo, all pairs are composed of satellite- 
central subhalos, while for Minf > 10 11 h~ 1 MQ, 25% of close 



• The radial profile of recently merged satellites roughly 
follows that of the entire satellite population, which also 
follows the halo density profile out to the virial radius. 
Satellite-satellite mergers preferentially occur in the outer 
regions of a halo. 
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• Mergers defined via major mass gain exhibit a strong 
small-scale clustering enhancement because significant mass 
gain is caused by interactions with a neighboring subhalo. 
However, in these cases the neighbor exists both before and 
after the 'merger' (i.e., no coalescence), and the major mass 
gain criterion does not correspond to dynamically disturbed 
subhalos. Thus, we do not expect these objects to correlate 
with active galaxies. 

Our measured subhalo merger enhancement suggests 
that, for populations with the same stellar mass at high 
redshifts, recent galaxy mergers should exhibit excess clus- 
tering at small radius, with a possible decrement in clus- 
tering at slightly larger scales. This effect will persist for 
only a short period of time after the merger, with a stronger 
excess signal for galaxies of higher stellar mass. Measuring 
this small-scale cross-correlation (against a general popula- 
tion) requires a surprisingly large volume (> 100/i -1 Mpc) 
because massive (satellite rich) halos contribute significant 
signal to the cross-correlation. Thus, a fair sample of rich 
halos is necessary for robust conclusions. 

A comparison of the clustering of our mergers to obser- 
vational phenomena is a nontrivial future step, since many 
merger observables depend on complex gas physics. Our re- 
sults do not include whether the subhalos are gas-rich or not, 
though at the high redshifts we examine, we expect almost 
all galaxies to be gas-rich. Sate llites can be stripped of much 
of th eir gas during infall (e.g.. iDolag et all feoOS: Sar o et all 
2008), though we have considered relatively massive satel- 
lites which have short infall times and thus experience less 
gas stripping. 

Time scales for observables after a merger also have 
large uncertainty and scatter. For example, an optical 
quasar might appear only ~ 1 Gyr afte r the merger (e.g., 
iHopkins et "aill2008l ; [Springel et alj|2005l ). by which time the 
enhanced clustering that we see has disappeared. X-ray sig- 
nals might appear sooner0 Starbursts or starburst-related 
objects (e.g., sub-millimeter galaxies, Lyman break galaxies, 
ULlR Gs) might also c ommence more rapidly after a merger 
(e.g., ICox et ail l2008h and thus appear a more promising 
analog to the effects we find here. 

While w e were preparing this work for preparation 
I Angulo ~et al. (2008) appeared which also considers the host 
halos of satellite mergers, in the Millennium simulation. 
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